clear all
mycell = cell(1, 4, 1)
mycell{1} = [0.5 0 1]'
mycell{2} = [-2 1 3]'
mycell{3} = [0 1.4 5]'
mycell{4} = [-1 0.5 2]'

for step = 1: 100,
cen = compute_grassmann_mean(mycell, step); 
v = [mycell{1,:}]

plot3(v(1,:), v(2,:), v(3,:), '.'); 
hold on; 
plot3(cen(1),cen(2),cen(3), 'x', 'color', 'red');
end